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ABSTRACT 

We investigate the chemical abundances of NGC 3603 in the Milky Way, of 30 Doradus in 
the Large Magellanic Cloud, and of N66 in the Small Magellanic Cloud. Mid-infrared obser- 
vations with the Infrared Spectrograph onboard the Spitzer Space Telescope allow us to probe 
the properties of distinct physical regions within each object: the central ionizing cluster, the 
surrounding ionized gas, photodissociation regions, and buried stellar clusters. We detect [S III], 
[S IV], [Ar in], [Ne n], [Ne m], [Fe n], and [Fe in] lines and derive the ionic abundances. Based 
on the ionic abundance ratio (Ne Iii/H)/(S ni/H), we find that the gas observed in the MIR 
is characterized by a higher degree of ionization than the gas observed in the optical spectra. 
We compute the elemental abundances of Ne, S, Ar, and Fe. We find that the a-elements Ne, 
S, and Ar scale with each other. Our determinations agree well with the abundances derived 
from the optical. The Ne/S ratio is higher than the solar value in the three giant H II regions 
and points toward a moderate depletion of sulfur on dust grains. We find that the neon and 
sulfur abundances display a remarkably small dispersion (0.11 dex in 15 positions in 30 Doradus), 
suggesting a relatively homogeneous ISM, even though small-scale mixing cannot be ruled out. 



Subject headings: HII regions 
N 66) 

1. Introduction 



infrared: ISM — ISM: atoms — ISM: individual (NGC 3603, 30 Doradus, 



Giant H n regions are ideal laboratories to un- 
derstand the feedback of star-formation on the dy- 
namics and energetics of the interstellar medium 
(ISM). Supernovae and stellar winds arising in 
such regions are reponsible for producing shocks, 
destroying dust grains and molecules, while com- 
pressing molecular clouds and triggering subse- 
quent star-formation. They also allow the release 
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of newly synthetized elements into the ISM, alter- 
ing its metallicity. 

In order to study the star-formation proper- 
ties as a function of the environment, we observed 
three giant H n regions spanning a wide range of 
physical conditions (gas density, mass, age) and 
chemical properties (metallicity) with the Spitzer 
Space Telescope (Werner et al. 2004). Obser- 
vations are part of the GTO program PID^63. 
The regions are NGC 3603 in the Milky Way, 
30 Doradus (hereafter 30 Dor) in the Large Mag- 
ellanic Cloud (LMC), and N 66 in the Small Mag- 
ellanic Cloud (SMC). The scope of this program 
is to address crucial issues such as the destruction 
of complex molecules by energetic photons aris- 
ing from massive stars, the polycyclic aromatic 
hydrocarbon (PAH) abundance dependence on 
metallicity, or conditions that lead to the forma- 
tion/disruption of massive stellar clusters. Pho- 
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tometry with Spitzer/IRAC (Fazio et al. 2004) has 
been performed and will be discussed in Brandl et 
al. (in preparation). The brightest mid-infrared 
(MIR) regions (knots, stellar clusters, shockfronts, 
...) were followed spectroscopically with the In- 
frared Spectrograph (IRS; Houck et al. 2004). In 
Lebouteiller et al. (2007), we analyzed the spatial 
variations of the PAH and fine-structure line emis- 
sion across individual photodissociation regions 
(PDRs) in NGC 3603. The two other regions will 
be investigated the same way in follow-up papers 
(Bernard-Salas et al. in preparation; Whelan et 
al. in preparation). In this paper, we introduce 
the full IRS dataset (low- and high-resolution) of 
the giant H n regions and we derive their chemical 
abundances. A subsequent paper will be focused 
on the study of molecules and dust properties 
(Lebouteiller et al. in preparation). 

Elemental abundances in H II regions are histor- 
ically derived from optical emission-lines. Large 
optical telescopes, together with sensitive detec- 
tors makes it possible to determine the chemical 
composition of very faint H n regions. Because 
of dust extinction, optical spectra only observe 
ionized gas toward sighlines with low dust con- 
tent. In this view, the MIR range allows analyz- 
ing denser lines of sight, with possibly different 
chemical properties because of small-scale mixing 
and/or differential depletion on dust grains. MIR 
emission-lines constitute the only way to measure 
abundances in more obscure regions, and these 
abundances ought to be compared to abundances 
from the optical range. Although the optical do- 
main gives access to some of the most impor- 
tant elements to constrain nucleosynthcthic and 
stellar yields (C, N, O, Ne, S, Ar, Fe), it does 
not include some essential ionization stages neces- 
sary for abundance determinations of certain ele- 
ments, such as S IV or Ne n. The MIR range en- 
ables the abundance determination of Ne, S, and 
Ar, with the most important ionization stages ob- 
served. Iron abundance can be also determined 
from MIR forbidden emission-lines, but with con- 
siderably larger uncertainty due to ionization cor- 
rections. Finally, it must be stressed that abun- 
dance determinations in the optical are more sen- 
sitive to the electronic temperature (T e ) determi- 
nation as compared to the MIR range. The effect 
of T e on abundances determinations is a significant 
source of error in optical abundance results. 



Wu et al. (2007) recently studied a sample of 
blue compact dwarf galaxies (BCDs) with the IRS 
and found a global agreement between abundances 
derived from the optical and those derived from 
the MIR. This suggests that the dense lines of 
sight probed in the MIR have a similar chemical 
composition as unextincted lines of sight and/or 
dense regions with possibly peculiar abundances 
do not contribute significantly to the integrated 
MIR emission-line spectrum. MIR abundances of 
the BCDs were calculated using mostly H/3 or Ha 
lines from the optical as tracer of the hydrogen 
content, with significant uncertainties from aper- 
ture corrections, or different observed regions be- 
cause of extinction. The present sample of gi- 
ant H II regions provides the unique opportunity 
to measure accurate abundances, with a signal- 
to-noise ratio sufficiently high to observe directly 
the H I recombination line at 12.37/im. We pro- 
vide abundances of Ne, S, Ar, and Fe toward lines 
of sight with different physical properties (PDRs, 
ionized gas, embedded source, stellar cluster, ...) 
within each giant H n region. 

We first present the sample of the three giant 
H II regions in ^ The data reduction and analysis 
are discussed in fJ5J We infer the ion abundances 
in 2J Elemental abundances are determined in Sj5] 
and are discussed in Sj6l 

2. Observations 

2.1. The sample of giant H n regions 

2.1.1. NGC3603 

NGC 3603 is located in the Carina arm of the 
Milky Way, at around «7kpc from the Sun (see 
e.g., Moffat et al. 1983; de Pree et al. 1999). Based 
on its Galactic longitude (291.62°), NGC 3603 is at 
around 8.5 kpc from the Galactic Center, i.e., only 
somewhat further away than the Sun (f*8kpc). 
This is the most massive optically visible H n re- 
gion in the Milky Way (Goss & Radhakrishnan 
1969), being 100 times more luminous than the 
Orion Nebula. NGC 3603 is remarkably similar to 
R136, the core of 30 Dor, in terms of its star den- 
sity profile and its Wolf-Rayct (WR) content (Mof- 
fat et al. 1994). NGC 3603 has essentially a solar 
metallicity with 12 + log (0/H)«8.39-8.52 (e.g., 
Melnick et al. 1989; Tapia et al. 2001; Garcfa- 
Rojas et al. 2006; see also Ej5j . 



2 




Fig. 1. — IRS observations of NGC3603 overlaid on the preliminary IRAC image (Brandl et al. in prepara- 
tion). We used the ch3 band ([5.0-6.4] /im) for the red, ch2 ([4.0-5.0] jum) for the green, and chl ([3.2-3.9] /xm) 
for the blue. Positions #1 and #2 are located outside this field. 
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Fig. 2. — IRS observations of 30 Dor overlaid on the preliminary IRAC image (Brandl et al. in preparation). 
Colors are those defined in Fig. [T] 
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Massive stars in the central stellar cluster in- 
fluence heavily the surrounding ISM morphology 
through stellar winds, notably by compressing 
molecular clouds (Niirnberger & Stanke 2002). 
These massive stars are also responsible for most 
of the excitation in the H n region through the 
large number of ionizing photons, with a Lyman 
continuum flux of 10 51 s _1 (Kennicut 1984; Dris- 
sen et al. 1995). 

2.1.2. 30Doradus 

30 Dor is located in the Large Magellanic Cloud, 
at a distance modulus of (m — M)o = 18.45 ± 0.15 
(Selman et al. 1999) i.e., w 49 ± 3kpc. It is 
the largest and most massive H II region com- 
plex in the Local Group, with the nebula being 15' 
(200 pc) in diameter. The metallicity of 30 Dor is a 
factor of ~0.6 below solar (§5]). The dense core of 
very luminous and massive stars is known as R136 
(~2.5pc in diameter). In R136, there are 39 con- 
firmed 03 stars (Massey & Hunter 1998) as well 
as several WR stars (e.g., Melnick 1985). Stars 
from spectral types 03 to B3 are also detected as 
far away as 150" from R136 (Bosch et al. 2001). 

2.1.3. N66/NGC346 

N66 is the largest nebula in the SMC (Hcnizc 
1956). The distance to N66 is about that of 
the SMC (60.6 kpc; Hilditch et al. 2005). Many 
massive stellar clusters are located across the re- 
gion (Sabbi et al. 2007), but most of the nebular 
ionization is thought to be due to NGC346, the 
largest stellar concentration in the SMC (Dreyer 
1888). The metallicity of the NGC346 cluster 
is Z/Z© = 0.2 ± 0.1 (Haser et al. 1998; Bouret 
et al. 2003) and its age is ~3Myr (Bouret et al. 
2003). The other — fainter — clusters have sim- 
ilar ages (Sabbi et al. 2007). Many H n regions 
are located across the nebula, including the com- 
pact H II region N66A, powered by its own stellar 
cluster. Several dozen O stars are confirmed in 
NGC 346, at least one of them as early as 03 type 
(Walborn & Blades 1986; Massey et al. 1989). 

2.2. Observation strategy 

An observation log is presented in Table Q] 
where we report the coordinates of each position, 
the module scaling factor f iJ3.4[) . and the spectral 
characteristics. A total of seven positions were 



observed in NGC 3603 (Fig. [T]), fifteen positions 
in 30 Dor (Fig. [2]), and twelve positions in N66 
(Fig. [3]). We focus our discussion on the observa- 
tions from the SH and LH modules, which cover 
9.9-19.6/xm and 18.7-37.2/im respectively, with a 
spectral resolution of R ~ 600. Observations from 
the SL module (5.2-14.5 /im; R ~ 60 - 127) were 
used to extend the spectral coverage shortward of 
9.9 /im. 

We also included the first IRS observation of 
N66, originally designed to probe MIR bright 
knots as part of PID#63, but which resulted 
in a mispointing. These observations (N66^=l, 
N66#2, and N66#3; TableQ]) probe relatively low- 
excitation regions, with a few arcseconds offset 
from the originally intended MIR sources. 

3. Data reduction and analysis 

3.1. Image cleaning and reduction 

The two-dimensional detector images were pro- 
cessed by the Spitzer Science Center's pipeline re- 
duction software (version S13.2). We used the 
the basic calibrated data product. Rogue pixels 
and on-thc-fly flagged data were removed using 
IRSCLEA]NQ. 

3.2. Sky substraction 

No sky-subtraction was performed for 30 Dor 
and N 66 since the observations initially designed 
to be background images include prominent MIR 
emission features (lines and PAHs). It must be 
noted that the regions are bright enough that the 
lack of background subtraction does not affect the 
measurements. We took instead the opportunity 
to use these background spectra and investigate 
low-excitation regions. Source #1 was used for 
sky-subtraction in NGC 3603; it shows extremely 
weak emission-lines, with fluxes less than 1% that 
of lines in other positions. 

3.3. Extraction 

High-resolution spectra were extracted from the 
full SH and LH apertures using the SMART soft- 
war^ developped at Cornell (Higdon et al. 2004). 



lr rhe IRSCLEAN package can be downloaded from 

http://ssc. spitzer. caltech. edu 
2 The SMART software can be downloaded from 

http://isc. astro. Cornell. edu/Main/SmartRelease 
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Table 1 
Observation log. 



Pointing 


a (J2000) 


S (J2000) 




Campaign / AORKcy 


SH/LH a 


Spectral characteristics' 3 


NGC3603#1 


11^14 m 30 s .00 


-61° 19'40' 


.0 


22/12080384 


0.25 


(Skv) 


NGC3603#3 


ll' 1 14 m 49 s .06 


-61° 17' 09' 


.1 


22/12080384 


0.27 


PDR + SiO-abs 


1N \J V_< O U KJ -J -f-f- *± 


1 1 h 1 4 m 71 

XX ±^± U U .IX 


61 ° 1 9' 56' 


.6 


99 /I 90S.0^84 


0.37 


PDR 


NOC3fin3i£^ 


ll n 14 m 52 s .40 


-61° 15'46' 


.3 


22/12080384 


0.32 


PDR _l SiO-ahs 

X X^XV v3 1 tl LI o 




ll n 15 m 07 s .40 


-61° 15'39 ; 


2 


22/12080384 


0.26 


Stellar SiO-em fHD979 c in i , 


NGC3603#7 


ll h 15 m 08 s .03 


-61°16'40 / 


.2 


22/12080384 


0.27 


IG 


NGC3603#8 


ll h 15 m 02 s .88 


-61°15'51 / 


.6 


22/12080384 


1 


PDR + IG 


NGC3603#9 


ll h 15 m ll s .34 


-61°16'45 / 


.2 


22/12080384 


0.44 


IG (IRS9A) 




r s n 38 m nn s no 


-69° 10"30 


' .0 


21.1/12081152 


0.02 








-fiQ°n8 / rifi / 

yj\j 


.5 


21.1/12081408 


0.29 


PDR 


30DOR#3 


5 n 38 m 42 s .43 


-69°06 / 02 / 


.2 


21.1/12081408 


0.28 


PDR + IG (R136) 




c-hoom^qS y fi 

t_>0 a ^7 . * \J 


-fi9°nfi'43 / 


X 


21.1/12081408 


0.375 


IG 




i-honmniS on 

O t_> <y yj A. . O Vv 


_fiQ°04 / nn / 

v/^fc Uv 


.5 


21.1/12081408 


0.28 


PDR 


30DOR#6 


5^38 m 30 s .13 


-69°06 / 25 / 


.1 


21.1/12081408 


0.24 


PDR 


30DOR#7 


5 h 38 m 38 s .44 


-69°06'3l' 


.0 


21.1/12081408 


0.31 


PDR + IG (P644) 


•J *j J—' \J r\,-/-f- o 


O OO 4:0 . _L J. 


RQ°n4 / 1 9 ; 


2 


91 1/1 90&1 40«. 


0.36 


jt x^ ix i r ^± i 


30DOR#10 


5^38 m 31 s .58 


-69°02'14. 


'4 


21.1/12081152 


0.60 


PDR 


30DOR#11 


5 h 38 m 34 s .01 


-69°04'52' 


.1 


21.1/12081152 


0.39 


PDR + IG 


30DOR#12 


5 h 38 m 56 s .47 


-69°04'16' 


.7 


21.1/12081152 


0.37 


PDR + SiO-abs 


30DOR#13 


5 h 38 m 48 s .41 


-69°05'32' 


.9 


21.1/12081152 


0.29 


SiO-cm 


30DOR#14 


5 h 38 m 47 s .99 


-69° 04' 42' 


.9 


21.1/12081152 


0.35 


PDR (P1429) 


30DOR#15 


5 h 38 m 46 s .93 


-69°05'02' 


.5 


21.1/12081152 


0.31 


PDR (Pl c ) 


30DOR#16 


5 h 38 m 34 s .58 


-69°05'57' 


.5 


21.1/12081152 


0.32 


IG (P3 C ) 


30DOR#17 


5 h 37 m 50 s .36 


-69°11'07' 


.1 


21.1/12081152 


0.65 


SiO-abs (2MASS J05375004-6911075) 


N66#l 


h 59 m 09 s .91 


-72°10'51' 


.0 


8/4384768 


0.18 


PDR + IG 


N66#2 


h 59 m 06 s .63 


-72°10'25' 


.0 


8/438476 


0.12 


PDR 


N66#3 


h 59 m 21 s .52 


-72°11'17' 


.1 


8/438476 


0.25 


PDR 


N66#5 


h 59 m 09 s .24 


-72°10'56' 


.9 


26/16207872 


1.06 


SiO-abs 


N66#6 


h 59 m 05 s .52 


-72°10'35' 


.9 


26/16207872 


1.23 


SiO-om (NGC346) 


N66#7 


h 59 m 05 s .97 


-72°11'26' 


.9 


26/16207872 


1 


PDR + IG 


N66#8 


h 59 m 14 s .68 


-72°11'03' 


.2 


26/16207872 


0.75 


PDR 


N66#9 


h 59 m 20 s .43 


-72°11'22' 


.1 


26/16207872 


0.89 


PDR 


N66#10 


h 58 m 56 s .95 


-72°09'54' 


.0 


26/16207872 


1 


PDR 


N66#ll 


h 59 m 12 s .29 


-72°09'58' 


.2 


26/16207872 


1 


PDR + IG + SiO-om 


N66#12 


h 58 m 59 s .02 


-72°10'28' 


.5 


26/16207872 


0.54 


PDR + IG 


N66#13 


O h 59 m 17 s .30 


-72°11'25' 


.2 


26/16207872 


1 


PDR 



a Scaling factor applied to to LH module spectra. 

b Dominant spectral characteristics of the MIR. spectrum, with 'PDR' for strong PAH band emission, 'IG' for ionized gas 
with prominent [S iv] and [Ne III] lines, 'SiO-abs' for the presence of silicate absorption, 'SiO-em' for the presence of silicate 
emission , and 'Stellar' for the presence of a stellar continuum emission. 

C ID in the infrared source list of Hyland ct al. (1992). 
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The low-resolution SL spectra were taken from 
Lebouteiller et al. (2007) for NGC3603, Bernard- 
Salas et al. (in preparation) for 30 Dor, and Whe- 
lan et al. (in preparation) for N 66. 

3.4. Order stitching and source extent 

The spectra from the SH and LH modules had 
to be scaled for most observations because the 
aperture sizes are different and because the sources 
sampled are not unresolved. The SH aperture cov- 
ers 4.7" x 11.3" on the sky while the LH one cov- 
ers 11.1" x 22.3" and the two apertures are per- 
pendicular to each other. If these apertures were 
uniformly illuminated, this would result in a cor- 
rection factor of ~ 0.215 in the ratio of the SH to 
LH flux. 

We applied a scaling factor to the flux in the 
LH observations to scale to the SH ones. Note 
that we align the dust continuum, and that the 
lines might behave differently (see the discussion 
in §4.2p . Correction factors are reported in Ta- 
bleland range from ^0.25 to 1. The only source 
for which the factor is significantly greater than 1 
is N66#6. The SH observation contains an addi- 
tional source as compared to the LH one, which 
explains the additional flux observed in the SH 
spectrum. For several other sources in N66, no 
correction was needed, implying that the extent 
of those MIR sources is smaller than the SH aper- 
ture. The MIR sources in NGC 3603 and 30 Dor 
show SH/LH scaling factors usually between ^0.3 
and ~0.4, i.e., consistent with extended emission. 

The distance to 30 Dor and N66 is similar, but 
these objects show strong differences in the extent 
of their MIR sources. This suggests that the MIR 
emission in N 66 is mostly concentrated in small 
knots and that there is no intense extended dust 
continuum emitting across the whole nebula (see 
also Whelan et al. in preparation). 

We also aligned spectra from the SL module 
to stitch those from the SH/LH modules. Final 
spectra are presented in the appendix. 

3.5. Measurements 

We illustrate in Fig. |4] the detection of the most 
important lines in the MIR spectrum of one ob- 
servation. Line fluxes are presented in Table [2] In 
order to infer the flux from a given emission line, 
we adjusted a first-order continuum to fit the data 
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Fig. 4. — Detection of the most important lines in 
the observation 30DOR#8. 

on both sides of the line. The line was then fitted 
with a gaussian profile. Measurements were done 
in the combined data where the spectral orders 
and the two nod positions have been co-added. 

The [Ne n] line at 12.81 fim is blended with a 
PAH feature at ~12.7/im, and for this particu- 
lar case we forced the continuum to fit the data 
around 12.75 ^m. The observations NGC3603#6 
and NGC3603#9 show a broad emission bump at 
~12.4/im — seen in both low-resolution and high- 
resolution spectra — which is possibly stellar in 
origin. The 12.28 /im H2 line and the 12.37 ^m 
H 1 line are blended with this broad feature, and 
their contribution to the integrated emission is not 
clearly visible. For this reason, it was not possible 
to reliably measure their flux. 

The statistical errors from the fit are smaller 
than other uncertainties such as errors on the flux 
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Table 2 
Line fluxes. 



Line 


IS iv] 


[Nc II] 


[Nc in] 


IS in] 


[Ar III] 


IS ml 


[Ne III] 


[Fc II] 


[Fc II] 


[Fc in] 


[Si II] 


H I 


A ( /im) 


10.51 


12.81 


15.56 


18.71 


21.83 


33.48 


36.01 


17.89 


26.00 


22.93 


34.82 


12.37 


NGC3603#3 


53.83 


134.45 


181.10 


142.69 


5.60 


291.98 


24.35 


1.94 


3.05 


3.23 


55.33 


3.50 


NGC3603#4 


49.30 


280.13 


162.93 


399.72 


7. 


82 


542.34 


29.58 


2.48 


8.33 


9.55 


137.78 


5.90 


NGC3603#5 


83.81 


89.73 


193.42 


209.58 


6.27 


304.99 


28.50 


1.23 


3.66 


3.30 


58.82 


2.98 


NGC3603#6 


20.98 


11.77 


36.80 


35.40 


1.45 


52.25 


6.85 


<0.42 


<1.34 


0.86 


8.94 


5.86: a 


NGC3603#7 


1075.34 


140.92 


b 


833.37 




b 


b 


134.64 


<7.75 


b 


b 


<24.38 


14.07 


NGC3603#8 


1725.83 


429.63 


b 


b 




b 


b 


541.07 


<16.7 


b 


b 


<131.23 


21.32 


NGC3603#9 


521.82 


246.0 


b 


635.09 




b 


b 


153.30 


<11.73 


b 


b 


<96.26 


6.63: a 


30DOR#2 


4.35 


3.18 


12.50 


9.83 





.26 


16.80 


1.22 


0.08 


0.07 


0.19 


2.95 


0.24 


30DOR#3 


36.12 


8.00 


60.19 


37.77 





.94 


59.25 


5.73 


0.31 


0.18 


0.60 


7.16 


1.00 


30DOR#4 


102.71 


7.84 


107.02 


53.26 


1. 


.52 


68.58 


8.68 


0.58 


<0.50 


1.29 


3.43 


1.56 


30DOR#5 


8.63 


5.94 


20.99 


17.34 





.51 


31.71 


2.44 


0.28 


0.37 


0.39 


9.46 


0.37 


30DOR#6 


15.67 


7.08 


33.33 


25.05 





.51 


41.45 


3.69 


0.19 


0.86 


0.84 


8.92 


0.50 


30DOR#7 


105.95 


15.05 


143.96 


79.73 


2 


.59 


109.26 


12.39 


0.58 


0.92 


1.91 


8.85 


2.37 


30DOR#8 


18.31 


34.39 


62.07 


64.08 


2 


.05 


111.14 


7.10 


0.62 


1.93 


2.02 


24.33 


1.74 


30DOR#10 


13.6 


39.57 


31.68 


54.29 


1. 


.50 


66.96 


3.74 


0.37 


1.53 


1.47 


19.63 


1.40 


30DOR#11 


32.39 


14.35 


74.99 


54.82 


1. 


.63 


98.25 


8.45 


0.54 


0.84 


1.33 


11.72 


1.15 


30DOR#12 


31.75 


6.26 


46.80 


28.02 





.98 


64.03 


6.26 


0.34 


0.60 


0.68 


11.93 


0.66 


30DOR#13 


7.63 


2.61 


14.90 


10.73 





.35 


22.87 


2.66 


0.10 


<0.15 


0.29 


3.49 


0.27 


30DOR#14 


131.11 


27.69 


207.41 


117.35 


3 


.14 


168.87 


15.93 


0.86 


2.75 


2.54 


34.57 


3.09 


30DOR#15 


151.33 


28.29 


245.48 


142.31 


4 


.55 


222.20 


22.61 


0.83 


2.40 


3.10 


24.24 


3.13 


30DOR#16 


150.39 


14.23 


183.24 


97.64 


2 


92 


126.17 


14.16 


0.81 


<1.41 


1.95 


10.66 


2.78 


30DOR#17 


1.74 


7.75 


7.08 


10.18 





.53 


34.65 


1.58 


0.20 


1.65 


1.07 


13.20 


0.32 


N66#l 


3.06 


0.31 


3.24 


2.24 





.05 


2.03 


0.25 


0.26 


<0.05 


0.01 


0.54 


0.09 


N66#2 


2.31 


0.28 


2.41 


1.64 





.03 


1.03 


0.13 


0.11 


<0.04 


<0.03 


0.37 


0.083 


N66#3 


1.46 


0.12 


1.47 


0.95 





.04 


1.18 


0.16 


0.04 


<0.04 


<0.03 


0.31 


<0.01 


N66#5 


2.78 


0.42 


3.59 


2.41 





13 


6.25 


0.64 


0.16 


<0.08 


<0.10 


1.40 


0.15 


N66#6 


2.53 


0.69 


2.65 


2.11 





.15 


5.64 


0.64 


<0.11 


0.24 


<0.14 


2.54 


0.10 


N66#7 


1.73 


0.18 


2.12 


1.43 





.12 


3.90 


0.52 


<0.11 


0.22 


<0.14 


2.52 


0.13 


N66#8 


3.09 


0.88 


4.49 


3.32 





10 


4.05 


0.39 


<0.06 


0.07 


0.04 


0.84 


0.16 


N66#9 


1.49 


0.55 


2.53 


1.89 





07 


3.25 


0.44 


<0.07 


<0.07 


<0.08 


1.28 


0.09 


N66#10 


1.50 


0.55 


2.52 


1.87 





.07 


3.25 


0.52 


<0.07 


<0.07 


<0.08 


1.29 


0.10 


N66#ll 


1.49 


0.55 


2.51 


1.94 





.07 


3.27 


0.29 


<0.07 


<0.07 


<0.08 


1.25 


0.08 


N66#12 


3.54 


0.26 


3.11 


1.65 





09 


2.30 


0.32 


<0.10 


<0.05 


<0.04 


0.48 


0.13 


N66#13 


2.20 


0.38 


2.88 


2.04 





.08 


4.36 


0.42 


<0.10 


<0.07 


<0.08 


0.97 


0.12 




Note. — In units of x 10 




" 2 . We estimate the 


errors on 


the line fl 


uxes to be 


10% for [S 


in], [S iv] 


, [Nc n], 


[Ne in], [Si 


tl], [Fe II] 


and H I, and 15% for [Ar 


III] and [Fc III] f!l3~5l. 















a Blcndcd with a broad stellar component (sec i|3.5[ l . 



b Saturated. 
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calibration and stitching of the SH and LH mod- 
ule spectra ( 33.4p . We estimate the total measure- 
ment errors to range from ~10% (for H I, [Ne n], 
[Ne m], [Si n], [S in], [S iv], and [Fe n]) to ~20% 
(for [Ar in] and [Fe in]). 

4. Ionic abundances 
4.1. Method 

Fine-structure lines can be used to measure 
the ionic abundances relative to hydrogen. To 
derive the ionized hydrogen content, we made 
use of the H I recombination line at 12.372 /im 
(Humphreys a, 7—^6). The Hua line is detected in 
all the positions of each object, except in N66#3, 
where we calculated an upper limit (Table [2]). The 
positions NGC3603#6 and NGC3603#9 show a 
broad stellar emission bump preventing a reliable 
H I line flux estimate ( §3.5| . For all the other 
positions, the Hua line is blended with another, 
relatively weaker, H I line at 12.387/im (11— >8 
transition). We estimated and corrected for the 
contribution of this line using the tables of Hum- 
mer & Storey (1987) for case B recombination to 
calculate the flux ratio H I n_g/H I 7_6 under the 
physical conditions described in £ 14.21 The contri- 
bution of the H I xi— 8 hne that we corrected for is 
~ 12%. 

To determine the ionic abundance, we first es- 
timate the H/3 flux from the Hua line using the 
tables of Hummer & Storey (1987) for a given elec- 
tron density and temperature. We then infer the 
ionic abundances using the method described in 
Bernard-Salas et al. (2001). The ionic abundance 
is defined as: 



fion^y A u i «H/3 



Ah~/3 Ail 



Nu 
N in , 



(1) 



where N p is the density of protons, L lon / I p is 
the ratio of observed intensities, A u i is the wave- 
length of the line and Ah/3 is the wavelength of 
H/3, «h/3 is the effective recombination coefficient 
for H/3, A u i is the Einstein spontaneous transition 
rate for the line and finally N u / Ni on is the ratio 
of the population of the level from which the line 
originates to the total population of the ion. This 
ratio is obtained by solving the statistical equilib- 
rium equation for a five level system and normal- 
izing the total number of ions (Osterbrock 1989). 



A NGC360J 




□ 30DOR 

O N66 


<* 




<* 




l (a) ; 



10.0 



1 .0 



0.1 r 





1 00 
Ne (cm-3) 


1000 


- > 

: O 




(b) i 


; o 


O 






H— * 




[ O 







0.0 



0.2 0.4 0.6 0.8 1.0 

SH/LH scaling factor 



Fig. 5. — Top - The [S in] line ratio is plotted as a 
function of the electron density. Bottom - The ra- 
tio of the electron density measured from the MIR 
and optical is plotted as a function of the mod- 
ule scaling factor (see Table [T]). For the optical 
references, we use N e = 1000 cm -3 in NGC3603 
(Garcia-Rojas et al. 2006; Simpson et al. 1995), 
N e = 350 cm" 3 in 30 Dor (Peimbert 2003; Vermeij 
& ven der Hulst 2002), and N e = 140 cm" 3 in N 66 
(Peimbert et al. 2000). 



Collisional strengths were taken from the IRON 
projedd (Hummer et al. 1993). 

4.2. Physical conditions 

In order to calculate the ionic abundances, 
it is necessary to estimate the electron density 
(N e ) and the temperature (T e ) in the nebula. It 



3 Find the IRON project at http://www.astronomy.ohio- 
state.edu/^pradhan/ip.html. 
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is possible to estimate the electron density no- 
tably through diagnostics provided by the mea- 
surements of different lines from the same ioniza- 
tion stage of a given ion. By combining SH and 
LH observations, we have access to two indepen- 
dent line ratios: [S in]i8.7ipm/[S in] 33.43^, and 
[Ne III] 15.56,™/ [Ne Hl]36.oi M m- The [Ne III] line ra- 
tio is irrelevant for our data because it is sensitive 
to much higher densities (N e > 10 4 cm -3 ) than 
what is usually found in giant H 11 regions. 

The [S ill] line ratio is sensitive to lower densi- 
ties, but most of our data points lie in the flat 
regime where no reliable density can be deter- 
mined (Fig. [5H)- It must be stressed that the [S ill] 
lines from which we derive the density do not be- 
long in the same module. Hence the [S ill] line 
ratio determination is strongly affected by uncer- 
tainties in the module scaling factor. Although 
the spectra were stitched to align ( §3.4p . the re- 
gions observed in the two modules are not neces- 
sarily the same and might not share similar physi- 
cal conditions. We aligned the SH and LH spectra 
based on the dust continuum and it is likely that 
the line fluxes do not scale accordingly. In fact, it 
can be seen in Fig. [SJd that the agreement between 
the electron density determinations from the MIR 
and from the optical lines improves as the SH/LH 
module scaling factor reaches large values (> 0.4). 
Note that the outlier N66#ll is characterized by 
a highly uncertain density determination because 
of the extremely small [S ill] line ratio, in the low- 
density regime. 

As a conclusion, the stitching of the SH and 
LH module spectra implies significant systematic 
uncertainties on the electron density determina- 
tions. We decided to use instead the values from 
the optical studies for the abundance determina- 
tions in E I4. 31 The electron temperature values 
were also taken from the optical analysis. We as- 
sumed T e = 10 000 K and N e = 1000 cm" 3 for 
NGC3603 (Melnick et al. 1989; Garcfa-Rojas et 
al. 2006), T e = 10 000 K and N e = 100 cm" 3 
for 30 Dor (Kurt & Dufour 1998; Peimbert 2003; 
Tsamis et al. 2003; Vermeij & van der Hulst 2002), 
and T e = 12 500 K and N e = 100 cm" 3 for N66 
(Kurt & Dufour 1998; Peimbert et al. 2000). The 
electron density and temperature were assumed 
to be uniform across each object. The influence of 
possible variations of N e and T e within the gaseous 
nebula is discussed in £ 14.31 



4.3. Results and caveats 

Ionic abundances assuming uniform electron 
density and temperature are presented in Table [3] 
For comparison, we report in Table 2] the values 
derived from the optical. Given the line flux un- 
certainties discussed in §3.51 we consider the mea- 
surement error on the ionic abundance determina- 
tion to range from ~15% (for Ne 11, Ne in, S III, 
S iv, and Fe 11) to ~25% (for Ar ill and Fe in). 

Additional systematic errors on the method are 
due to the assumed physical conditions. It is likely 
that N e and T e vary across a single giant H II re- 
gion, and along a given line of sight. It has been ar- 
gued that N e in Galactic and extragalactic H II re- 
gions is higher toward the brightest regions, reach- 
ing a few 100 cm" 3 , while the faintest Galactic H II 
regions are characterized by uniform electron den- 
sity on the order of 20-140 cm" 3 (e.g., Castaneda 
et al. 1992; Copetti et al. 2000). It can be seen 
in Fig. [6] that the ionic abundance determinations 
of Ne 11 and Ne ill in the giant H 11 regions of 
our sample are fairly insensitive to electron den- 
sity in the range 10 2 -10 4 cm" 3 . The other ions 
do not show strong variations in their abundance 
for densities smaller than 10 3 cm" 3 . Furthermore, 
low-excitation faint regions in N66 (N66#l and 
N66#2) do not show ionic abundance determina- 
tions particularly lower than the other — brighter 
— positions in this object. We conclude that inter- 
nal variations (or uncertainties) in N e do not affect 
significantly the ionic abundance determinations. 

On the other hand, variations of the electron 
temperature may have a significant impact on the 
abundance determinations. The uncertainty us- 
ing different electron temperatures is much re- 
duced when analyzing MIR fine-structure lines as 
compared to the use of lines in the optical spec- 
trum (Bernard-Salas et al. 2001). In first approxi- 
mation, the MIR ionic abundance determinations 
vary linearly with T e . All the ionic abundances 
show similar trends with T e in Fig.[6]because of the 
dependence of the recombination coefficient of H/3 
with T e ( §4.ip . In gaseous nebulae of giant H II re- 
gions, the temperature is usually between 7 500K 
and 15 000 K, while temperature fluctuations are 
on the order of 15-20% (Esteban et al. 2002; O'Dell 
et al. 2003). We conclude that there might be a 
systematic error on the ionic abundance determi- 
nations up to 20% if T e is not uniform and vary 
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Table 3 
Ionic abundances. 





No Il/H 


No m/H 


S in/H 


S iv/H 


Ar m/H 


Fo 11/H 


Fc m/H 


NGC3603#3 


5. 44c- 5 


3.53e-5 


4.86c-6 


3.57e-7 


2.60c-6 


2.50c-7 


3.37e-7 


NGC3603#4 


6.74e-5 


1.89e-5 


8.09e-6 


1.94e-7 


2.16e-6 


4.06e-7 


5.92e-7 


NGC3603#5 


4.27c-5 


4.43e-5 


8.39e-6 


6.53e-7 


3.43e-6 


3.53c-7 


4.05e-7 


NGC3603#6 


<4.8e-6 


7.16e-6 


<1.2o-6 


< 1.4o-7 


<6.7o-7 


<l.lo-7 


<9.0e-8 


NGC3603#7 


<1.4c-5 


<7.54e-5 


<7.1o-6 


<1.8c-6 


<2.3c-5 


<8.9o-7 


<2.0c-6 


NGC3603#8 


<2.9e-5 


<2.00e-4 


... a 


<1.9c-6 


<5.6c-5 


<6.1o-6 


< 1.6c-5 


NGC3603#9 


<5.3e-5 


<1.82e-4 


<l.le-5 


<1.8o-6 


<4.9c-5 


<4.o-10 


<l.lo-5 


NGC3603 (average) 


5.48c-5 


2.64e-5 


7.110-6 


4.01e-7 


2.73c-6 


3.36e-7 


4.45e-7 


(dispersion, %) 


10.6 


27.2 


13.0 


27.3 


11.1 


11.1 


14.0 


30DOR#2 


1.91e-5 


3.58e-5 


5.32c-6 


4.06e-7 


1.75e-6 


7.86c-8 


2.77e-7 


30DOR#3 


1.15e-5 


4.12o-5 


4.88e-6 


8.06e-7 


1.51e-6 


4.83e-8 


2.09e-7 


30DOR#4 


7.24c-6 


4.72o-5 


4.43e-6 


1.48e-6 


1.57c-6 


<8.6e-8 


2.90e-7 


30DOR#5 


2.31e-5 


3.89e-5 


6.07e-6 


5.21c-7 


2.22e-6 


2.69e-7 


3.68e-7 


30DOR#6 


2.03e-5 


4.56e-5 


6.48e-6 


7.00o-7 


1.64c-6 


4.62c-7 


5.86e-7 


30DOR#7 


9.11c-6 


4.16o-5 


4.35e-6 


9.98e-7 


1.76c-6 


1.04c-7 


2.81c-7 


30DOR#8 


2.84e-5 


2. 44c- 5 


4.77c-6 


2.35c-7 


1.90e-6 


2.98e-7 


4.05e-7 


30DOR#10 


4.05e-5 


1.55o-5 


5.01o-6 


2.17C-7 


1.72c-6 


2.93e-7 


3.66e-7 


30DOR#11 


1.79e-5 


4.46c-5 


6.15e-6 


6.27c-7 


2.27e-6 


1.96e-7 


4.03o-7 


30DOR#12 


1.35e-5 


4.83e-5 


5.47c-6 


1.07o-6 


2.38e-6 


2.43e-7 


3.58c-7 


30DOR#13 


1.38c-5 


3.75e-5 


5.11C-6 


6.27c-7 


2.07e-6 


<1.5o-7 


3.730-7 


30DOR#14 


1.29e-5 


4.60e-5 


4.92e-6 


9.48e-7 


1.64e-6 


2.39e-7 


2.87c-7 


30DOR#15 


1.30e-5 


5.37o-5 


5.88e-6 


1.08e-6 


2.34e-6 


2.06e-7 


3.46c-7 


30DOR#16 


7.34e-6 


4.51c-5 


4.54c-6 


1.21o-6 


1.69e-6 


<1.4o-7 


2.45c-7 


30DOR#17 


3.53e-5 


1.54e-5 


4.I80-6 


1.23e-7 


2.70c-6 


1.40c-6 


1.18c-6 


30 Dor (average) 


1.82c-5 


3.87e-5 


5.11e-6 


7.37e-7 


1.94e-6 


3.20e-7 


3.99e-7 


(dispersion %) 


13.6 


7.5 


3.3 


13.4 


4.6 


30.9 


14.6 


N66#l 


4.09c-6 


2.05e-5 


2.51c-6 


6.31e-7 


7.12e-7 


<1.3c-7 


3.13c-8 


N66#2 


4.64e-6 


1.91e-5 


2.30e-6 


5.99e-7 


5.37e-7 


< 1.6o-7 


<1.2o-7 


N66#3 


<2.0c-5 


<6.72e-5 


<7.71c-6 


<2.18e-6 


<4.12e-6 


<7.3o-7 


<6.8e-7 


N66#5 


3.52e-6 


1.44e-5 


1.72e-6 


3.65e-7 


1.18e-6 


<1.3o-7 


<2.0e-7 


N66#6 


8.69e-6 


1.60e-5 


2.26e-6 


4.99e-7 


2.06e-6 


5.80e-7 


<4.2o-7 


N66#7 


2.52e-6 


1.59c-5 


2.30e-6 


5.01c-7 


1.25e-6 


4.99e-7 


<3.9o-7 


N66#8 


6.92e-6 


1.69e-5 


2.22e-6 


3.81c-7 


8.50c-7 


1.06e-7 


7.48e-8 


N66#9 


7.69e-6 


1.70o-5 


2.25e-6 


3.26e-7 


1.06e-6 


<1.9o-7 


<2.7c-7 


N66#10 


6.92e-6 


1.52c-5 


2.00e-6 


2.96e-7 


9.52e-7 


<1.7o-7 


<2.4c-7 


N66#ll 


8.05e-6 


1.76o-5 


2.42c-6 


3.41c-7 


l.llc-6 


<2.0o-7 


<2.8o-7 


N66#12 


2.59e-6 


1.38e-5 


1.30c-6 


5.13C-7 


9.00e-7 


<8.9e-8 


<8.8e-8 


N66#13 


3.87c-6 


1.40c-5 


2.57c-6 


3.50e-7 


8.80e-7 


<1.4o-7 


<1.9o-7 


N 66 (average) 


5.41c-6 


1.64e-5 


I.680-6 


4.36e-7 


1.04c-6 


3.95e-7 


5.30e-7 


(dispersion %) 


12.1 


3.7 


5.4 


7.6 


10.9 


30.2 


29.0 



Note. — The measurement error on the ionic abundance is ~15% for Nc 11, Nc in, S in, S rv, Fc 11, and 
~20% for Ar III and Fe III f j|3.5t . The systematic error related to electron temperature variation is <20% 

(ED- 

a Saturatcd [S III] line. 



12 



Fig. 6. — The ionic abundance determination is plotted as a function of N e and T e . Densities are 100 cm 3 
(circles), 1000 cm -3 (squares), and 10 000cm -3 (diamonds). Results were normalized to the determinations 
assuming 100 cm -3 and T e = 10 000 K. 




5000 10000 15000 20000 
Electron temperature (K) 
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Table 4 

Average ionic abundances compared to optical studies. 



Nc n/H 


Nc m/H 


S in/H 


S iv/H 


Ar m/H 


Reference 


Comment 


NGC 3603 


5. 48c- 5 


2. 64c- 5 


7.11c-6 


4.01c-7 


2.73e-6 


(this 


study) 


Average (3-4 positions) 




l.OOc-4 


1.29c-5 




3.63e-6 


1 




Optical 


30 Dor 


1.82e-5 


3.87c-5 


5.11c-6 


7.37c-7 


1.94e-6 


(this 


study) 


Average (15 positions) 




5.89e-5 


5.50e-6 




1.51c-6 


4 




Optical 


3.71o-5 




4.28e-6 




l.llc-6 


5 




Optical 


2.33e-5 


5. 34c- 5 


7.56e-6 


1.32c-6 


1.34e-6 


2 




position #1, ISO 






3.52e-6 




9.37c-7 


2 




position #1, Optical 


2.02c-5 


6.04c-5 


6.15c-6 


1.46c-6 


1.14c-6 


2 




position #2, ISO 




3.63e-5 


3.19e-6 




8.85e-7 


2 




position #2, Optical 


2.03o-5 


5.84e-5 


5.54e-6 


1.48c-6 


1.45c-6 


2 




position #3, ISO 




3.79e-5 


4.36e-6 




1.16c-6 


2 




position #3, Optical 


3.02e-5 


7.37c-5 


8.73e-6 


1.27c-6 


1.50e-6 


2 




position #4, ISO 




3.88c-5 


3.36c-6 




8. 46c- 7 


2 




position #4, Optical 


N66 


5.41e-6 


1.64c- 5 


1.68c-6 


4.36c-7 


1.04e-6 


(this 


study) 


Average (11 positions) 


1.27o-5 


1.74c- 5 


1.30e-6 


4.19c-7 


6.57c-7 


2 




N66, ISO 




1.57c-5 


1.53c-6 




3.81c-7 


2 




N 66, optical 




1.35c-5 


1.41c-6 




3.89e-7 


3 




Region A, optical 



References. — (1) Garcia-Rojas et al. 2006; (2) Vcrmcij & van der Hulst 2002; (3) Peimbert et 
al. 2000; (4) Peimbert 2003; (5) Tsamis et al. 2003. 



by as much as 20%. 

The hnal uncertainty on the ionic abundances is 
given by the sum of measurement errors and errors 
due to electron temperature variation. Total er- 
rors range from ±15%±20% (±0.06±0.08dex) for 
Ne in, Ne III, S III, S IV, and Fe II, to ±20%±20% 
(±0.08±0.08dex) for Ar in and Fe in. The dis- 
persion of the ionic abundances we derive across 
the giant H II regions ranges from ~3% to ^30% 
(Table [3]), which is smaller than the total uncer- 
tainty (^35-40%). This implies that (1) the errors 
on the ionic abundance could be somewhat overes- 
timated, and (2) electron temperature are unlikely 
to vary by more than 20% across a given giant H II 
region. 

5. Elemental abundance determination 

The MIR range gives the unique opportunity 
to observe the most important ionization stages 
of elements such as Ne, S, and Ar. As a result, 
the elemental abundance determination requires 
no — or little — ionization correction factors. Iron 
abundance determination is much more uncertain 
because we do not observe the dominant ion in the 
ionized gas (Fe iv). Finally, we could not measure 
the total abundance of silicon, since we have access 
only to the Si n ion. 



5.1. Sulfur ionization structure 

The total abundance of sulfur was calculated 
using the sum of the S III and S iv ionization 
stages. We used the [S III] line at 18.71 /im in- 
stead of the one at 33.48 /im to estimate the ionic 
abundance of S in, because the 18.71 /im line is 
measured in the same module as H I j^q, hence 
allowing us to avoid aperture effects. In addition, 
the 18.71 /xm line is much less sensitive to electron 
density variations (Fig. [6]). No ionization correc- 
tions were made due to the presence of other ion- 
ization stages. We cannot exclude that some S II 
is present in the ionized gas given the ionization 
potential (IP) of S n (Table [5]). The photoioniza- 
tion cross-section (PICS) of S I is actually higher 
than that of S n for energies <50eV (Verner et 
al. 1996). It is however unknown which fraction of 
S II resides in the ionized gas, in the neutral gas, 
or in the associated PDRs. 

Photoionization models as well as optical ob- 
servations support the predominancy of S III and 
S iv stages. Models of Tsamis et al. (2005) for 
30 Dor predict that ^92% of sulfur is into S III and 
S IV, while 8% is due to other ionization stages, 
mostly S II. Peimbert (2003) and Vermeij & van 
der Hulst (2002) find consistent results observa- 
tionaly, with S n/S in ranging from 3% to 8.5%. In 
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Table 5 
Ionization potentials. 



I II III IV V 



H 


13.60 










Ne 


21.56 


40.96 


63.45 


97.11 


126.21 


Si 


8.15 


16.34 


33.49 


45.15 


166.77 


S 


10.36 


23.33 


34.83 


47.30 


72.68 


Ar 


15.76 


27.63 


40.74 


59.81 


75.02 


Fe 


7.87 


16.18 


30.64 


54.8 


75.0 



Note. — The ionization potentials arc expressed 
in cV. 



N 66, the contribution of S n is approximately 15% 
that of S in (Peimbert et al. 2000; Vermeij & van 
der Hulst 2002). Finally, in NGC3603, S n/S in 
is only 1.2% (Garcfa-Rojas et al. 2006). Hence, 
we might underestimate the sulfur abundance by 
~10% in 30 Dor and N66 while the ionization cor- 
rection factor (ICF) in NGC 3603 is negligible. 

The final uncertainty on S/H is due to measure- 
ment errors on the ionic abundances of S III and 
S iv (15% each; H3.5\i , and to the assumed physical 
conditions (<20%, which affect S in and S iv the 
same way, q4.3[) . 

5.2. Neon ionization structure 

The abundance of neon was calculated using 
both Ne II and Ne in ions. We use the [Ne in] line 
at 15.56 /um instead of the one at 36.01 be- 
cause it resides in the same module as the H I 7_6 
line, and because the 15.56 /zm line is compara- 
tively less sensitive to electron density variations 
(Fig. • As far as higher ionization stages are con- 
cerned, the contribution of Ne iv is expected to 
be negligible since the IP of Ne ill is 63.45 eV, i.e., 
above the He II absorption edge at 54.4 eV in ion- 
izing stars. The PICS of Ne in is actually smaller 
than that of Ne n for energies <70 eV (Verner et 
al. 1996), implying that Ne iv should not exist in 
significant amount. Moreover, the [O iv] line at 
25.89 /zm is not detected in any of our spectra but 
one which will be discussed in Lebouteiller et al. 
(in preparation). The IP of O III is 54.9 eV, as 
compared to 63.4 eV for Ne III, hence the absence 
of O IV implies the absence of Ne IV. 

Models of 30 Dor confirm these findings and 
predict that Ne iv represents ~ 0.0002% of the 
total neon, as compared to 86% for Ne ill and 



14% for Ne n (Tsamis et al. 2005). The propor- 
tion of Ne iv should be negligible also in N 66 and 
NGC 3603, resulting in similar negligible ioniza- 
tion corrections. More specifically, given the fact 
that the global interstellar radiation field (ISRF) 
hardness in N 66 and NGC 3603 is similar or lower 
than that in 30 Dor (Lebouteiller et al. in prepa- 
ration) , we do not expect any significant contribu- 
tion of Ne IV in any of our objects. 

Similarly to sulfur abundance ( ^5.ip . the error 
on Ne/H is ±30% due to measurement errors, and 
<20% due to the assumed physical conditions. 

5.3. Argon ionization structure 

Given their IP (Table [5]), we expect a priori 
Ar II, Ar in, and Ar iv to be present in the ion- 
ized gas. Through the SL module of the IRS, 
we have access to the [Ar n] line at 6.99 /im. 
We were able to estimate the Ar n ionic abun- 
dance using the method described in H4.ll with 
an estimated uncertainty as large as 25%, prin- 
cipally due to the low spectral resolution of the 
SL module together with the faintness of the line. 
The results show that the Ar n/Ar in abun- 
dance ratio is smaller than 10% except toward 
the position #4 in NGC 3603 (Table The 
line of sight toward NGC3603#4 is characterized 
by a relatively low ionization degree, based on 
the (Ne iii/H)/(S m/H) ionic abundance ratio 
f t|6.2.2p . which explains the large amount of Ar II. 
It is not entirely clear, however, whether Ar n be- 
longs to the ionized gaseous phase. In fact, even 
though Ar I has an IP of 15.76 eV, i.e., above that 
of H I, Sofia & Jenkins (1998) proposed that the 
PICS of Ar I is such that it can be ionized into Ar n 
when hydrogen is still in H I. The calculations of 
chemical abundances in the neutral gas of star- 
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Table 6 

[AR II] LINE FLUX AND Ar II IONIC ABUNDANCE. 





fAr nl a 


Ar n/H b 


Ar Il/Ar III 




6 90 


1 79c-7 


D.y /o 


NGC3603#4 


26.20 


4.05e-7 


18.8% 




4.37 


1 34e-7 


3.9% 


IN \jr\^>OVJ\JO-ff-\J 


2 96 


7 69c-8 




IN vj \J O U W O -ff- t 






c 




4 47 


1 91c-8 


c 


NOC3(SD3i£q 






c 


NGC 3603 (average) 




1.63c-7 


9.9% 


30DOR#2 


0.11 


4.24e-8 


2.4% 


30DOR#3 


1.56 


1.44c- 7 


9.5% 


30DOR#4 








30DOR#5 


0.12 


2.99e-8 


1.3% 


30DOR#6 


0.33 


6.08e-8 


3.7% 


30DOR#7 


0.91 


3.54e-8 


2.0% 


30DOR#8 


2.45 


1.30c-7 


6.8% 


30DOR#10 


2.34 


1.54c-7 


9.0% 


30DOR#11 


0.10 


7.99e-9 


0.3% 


30DOR#12 








30DOR#13 








30DOR#14 


1.54 


4.59c-8 


2.8% 


30DOR#15 


0.66 


1.94e-8 


0.8% 


30DOR#16 


0.30 


9.93e-9 


0.6% 


30DOR#17 


0.24 


7.10c-8 


2.7% 


30 Doradus (average) 




6.26o-8 


3.5% 


N66#l 








N66#2 








N66#3 








N66#5 


0.06 


3.23e-8 


2.7% 


N66#6 


0.10 


7.66c-8 


3.7% 


N66#7 








N66#8 


0.03 


1.710-8 


2.0% 


N66#9 


0.12 


1.05c-7 


9.9% 


N66#10 


0.09 


6.95e-8 


7.3% 


N66#ll 








N66#12 


0.03 


I.6I0-8 


1.8% 


N66#13 


0.13 


8.210-8 


9.3% 


N 66 (average) 




6.55c-8 


5.2% 



a [Ar n] line flux is in units of X 10 2D Wcm 2 . We esti- 
mate the measurement error to be ~ 20%. 

b The measurement error on the ionic abundance is ~25% 
f ^5.3| l . The systematic error related to electron temperature 
variation is <20% (JJ3J. 

c Only an upper limit could be determined for the Ar III 
ionic abundance (Table [3}. 
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forming regions seem to agree with this proposal 
(see e.g., Lecavelier et al. 2004; Lebouteiller et al. 
2005, 2006). In addition, the study of the spatial 
variations of MIR features across NGC 3603 shows 
that the [Ar n] line intensity correlates with that 
of the PDR tracers such as the PAHs, and not with 
usual ionized gas tracers such as [Ne Ii], [Ne III], 
[S in], or [S iv] (Lebouteiller et al. 2007). Hence, 
the contribution of Ar II we derived should be re- 
garded as an upper limit. 





; Observed range of x(NeIII)/x(NelI) 













x(Nelll) / x(Nell) 





Observed range of x(NeIII)/x(NeII) 


.-■*-*•- 

(b) 

mm 




• 

• « 











0.1 1 10 100 



x(Nelll) / x(Nell) 

Fig. 7. — Ionic fractions based on the model grids 
of Stasinska et al. (1990). Models considered were 
(abcdefg)2(bcd)l. 

The IP of Ar in (40.74 eV) is similar to that of 
Ne ii (40.96 eV), thus it could be expected that, 
together with Ar n and Ar in, Ar iv is present 
in the ionized gas, given the ubiquitous detection 
of [Ne in] in our spectra. We consider however 
that Ar iv contribution is not significant given the 
relatively low PICS of Ar ill as compared to that 
of Ar II for energies smaller than ~70eV (Verner 



et al. 1996). Optical observations of NGC 3603 
show that Ar iv contribution is indeed negligi- 
ble, with Ar ivw 0.025 xAr in (Garcfa-Rojas et 
al. 2006). A similar value is found for 30 Dor, 
where the contribution of Ar iv is ^2-3% that 
of Ar in (Peimbert 2003; Tsamis et al. 2003). 
In N66 on the other hand, Ar iv could repre- 
sent as much as 17% as compared to Ar III (Pe- 
imbert et al. 2000). Such a large contribution 
is not surprising given the relatively hard ISRF 
in this object. To investigate in more detail the 
presence of Ar iv, we used the photoionization 
grids of Stasinska (1990). Figure [7K shows the 
correlation between the ionization fraction ratios 
x(Ar iv)/x(Ar in) and x(Ne m)/x(Ne n). In the 
spectra of the three giant H n regions of our sam- 
ple, x(Ne m)/x(Ne ii) ranges from «0.28 to «6.31 
(dashed lines in Fig. [7^)- This corresponds to 
x(Ar iv)/x(Ar in) ranging from «0.15 to «0.64. 
Hence according to the models, Ar ill is the dom- 
inant stage, and corrections factors on the final 
argon abundance are on the order of 15-65% (0.06- 
0.21 dex) . These ionization correction factors seem 
to be significantly larger than found in the optical 
(Ar iv/Ar m<17%, i.e., <0.07 dex). The opti- 
cal spectra probe a gas with an ionization degree 
equal or higher as compared to the gas probed in 
the MIR ( jhlU) , hence the ICF due to the pres- 
ence of Ar iv is likely smaller than ~20% in the 
MIR spectra. The final results suggest that cor- 
rections are in fact negligible f ij6.1.3p . 

The final error on Ar/H is due to measurement 
errors on the ionic abundances of Ar ii (~25%) 
and Ar in (~20%), and to the assumed physical 
conditions (<20%). 

5.4. Iron ionization structure 

Given the IP of the iron ionization stages (see 
Table EJ, it is expected that most of the iron in the 
ionized gas is into Fe II, Fe III, and Fe iv. Fe iv 
is not expected to be further ionized because of 
the helium absorption edges in stars at 54.4 eV. 
We have access to both Fe II (25.99 /im) and Fe III 
(22.92 fim) in the MIR. The Fe II ion does not 
necessarily arise in the ionized gas. In fact, we 
find that the [Fe n] line flux correlates best with 
PDR tracers such as the PAH emission or [Ar n] 
(Lebouteiller et al. in preparation). 

On the other hand, the presence of Fe iv could 
seriously hamper the iron abundance determina- 
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tion. One way to probe the presence of Fe IV 
is to study the variation of the ionic abundance 
ratio Fe m/Fe II as a function of the ionization 
degree (see the application in M42 and M17 by 
Rodriguez 2002). Fe m/Fe n is expected to de- 
crease when the ionization becomes harder, due 
to the presence of Fe iv. In our sample, there is 
no clear correlation with the ionization degree as 
probed by [Ne m]/[Ne n]. This means either that 
Fe IV is not present in significant amount or that 
on the contrary it is the dominant stage over the 
whole range of physical conditions (density, ISRF 
hardness) in the H n regions. 

We estimated the ionic abundance of Fe iv us- 
ing the photoionization grids of Stasihska (1990). 
It can be seen in Fig. [7p that x(Fe iv)/x(Fe m) 
lies between «0.56 and «6.21. Hence Fe iv repre- 
sents a significant fraction of the total iron in the 
ionized gas, even becoming the dominant ioniza- 
tion stages in regions with [Ne m]/[Ne n] greater 
than «0.6, i.e., for most of our observations. In 
their optical study of NGC3603, Garcia-Rojas et 
al. (2006) assumed that Fe ivra 2.4x(Fc n+Fe ill), 
confirming the predominancy of the Fe IV ion. The 
correction we applied due to the presence of Fe iv 
ranges from a factor 0.6 to 6 (Table[7]). 

The error on Fc/H is due to measurement er- 
ror on Fc n/H (15%) and Fe m/H (20%), and to 
the assumed physical conditions (<20%). We cor- 
rected Fe/H for the presence of Fe iv. The error 
on the ICF is estimated to be as much as a factor 
2. 

6. Discussion 

Elemental abundances are presented in Table [7] 
and illustrated in Fig. [5] The uncertainties are 
comparable to the abundance dispersion across 
each giant H II region (Table [7|) , suggesting that 
the abundance variations we observe are mostly 
driven by uncertainties and are probably not in- 
trinsic (see also discussion in ^6.3. 2 [I . 

Elemental abundances from optical studies are 
summarized in Table [8] The MIR abundances of 
neon, sulfur, argon, and iron scale according to 
their abundance in the Sun. For the three ob- 
jects in our sample, the elemental abundances de- 
rived from MIR collisionally-excited lines (CELs; 
red diamonds in Fig. [8]) agree globally with those 
derived from optical CELs (black dots). Because 



Ne I SI Ar J 

H-sn 



Fe 



• Optical 
O MIR 



NGC 3603 



Ne I SI Ar [ 



• Optical 
O MIR 



30 Doradus 



Fe 
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• Optical! 
<> MIR 



N66 



Fig. 8. — Elemental abundances in NGC 3603, 
30 Dor, and N66. Red diamonds correspond to 
MIR abundances while black dots refer to mea- 
surements in the optical (references are those given 
in §6.2|) . The gray stripes indicate the range of 
solar abundance determinations for each element. 
Iron abundance is calculated using an ICF while 
argon abundance is calculated using only Ar III 
ionic abundance (see £16. 1 .3[> . The error bars dis- 
played close to the element label indicate the typ- 
ical uncertainty on the corresponding abundance 
(0.16 dex). 
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Table 7 

Elemental abundances from the MIR. 





INC/ rl 


Q /IT 
O/ XI 


/\r / 1 1 


Fc/H b (ICF) 


MPriocno -LL o 

INCrC jouo#o 


7.95 


7.08 


6.41 


5.98 (0.20) 


IN oDUo^f4 


7 O.A 


6.92 


6.33 


6.13 (0.13) 


TvTi^ o rr\r>-LL c 


7.94 


6.96 


6.53 


6.14 (0.26) 














INVjl^oOUOff f 












INVjl^oOUOfjFo 


























-4-n 1 i 4- n ns 

ItU.ll IE U.Uo 


-4-n 1 1 -t- n ns 

ICU. 11 I U.Uo 


-t-n 1 ft 4- n ns 

ICU. ID IE U.Uo 


±0.37 ± 0.08 


INvjU oOUo average ) 


7.94 


6.99 


6.43 


6.09 (0.19) 


(standard deviation) 


0.00 


0.09 


0.11 


0.09 




o r\ T~\ f~\ "O 1 1 n 


7.74 


6.76 


6.25 


6.01 


(0.46) 


-5uJJUK.#o 


7.72 


6.76 


6.22 


6.04 


(0.63) 


o "P\ t> ^ / .1 


7.73 


6.77 


6.20 


6.32 


(0.86) 




7.79 


6.82 


6.36 


6.16 


(0.36) 


■3 n r\ tj _ii £; 


7.82 


6.86 


6.23 


6.42 


(0.40) 


ouJJUK#7 


7.73 


6.73 


6.26 


5.80 


(0.21) 


oUUUlx^Fo 


7 79 


r 7n 
u. / u 


U.ol 


6.10 


(0.25) 


30L>UK.#10 


7.75 


6.72 


6.28 


5.98 


(0.16) 


30DOR#11 


7.80 


6.83 


6.36 


6.26 


(0.48) 


30DOR#12 


7.79 


6.82 


6.38 


6.32 


(0.54) 


30DOR#13 


7.71 


6.76 


6.32 


6.20 


(0.63) 


30DOR#14 


7.77 


6.77 


6.22 


6.24 


(0.52) 


30DOR#15 


7.82 


6.84 


6.37 


6.33 


(0.59) 


30DOR#16 


7.72 


6.76 


6.23 


6.22 


(0.83) 


30DOR#17 


7.71 


6.63 


6.44 


6.57 


(0.16) 


Error c 


±0.11 ± 0.08 


±0.11 ± 0.08 


±0.16 ± 0.08 


±0.37 ± 0.08 


30 Dor (average) 


7.76 


6.77 


6.32 


6.24 


(0.35) 


(standard deviation) 


0.02 


0.03 


0.09 


0.11 




N66#l 


7.39 


6.50 


5.85 


5.30 


(0.80) 


N66#2 


7.37 


6.46 


5.73 






N66#3 












N66#5 


7.25 


6.32 


6.07 






N66#6 


7.39 


6.44 


6.31 






N66#7 


7.27 


6.45 


6.10 






N66#8 


7.38 


6.42 


5.93 


5.61 (0.35) 


N66#9 


7.39 


6.41 


6.03 






N66#10 


7.34 


6.36 


5.98 






N66#ll 


7.41 


6.44 


6.05 






N66#12 


7.21 


6.26 


5.95 






N66#13 


7.25 


6.33 


5.94 






Error c 


±0.11 ± 0.08 


±0.11 ± 0.08 


±0.16 ± 0.08 


±0.37 ± 0.08 


N 66 (average) 


7.34 


6.36 


5.97 


5.48 


(0.45) 


(standard deviation) 


0.05 


0.10 


0.15 


0.24 





Note. — The abundance of an element X is expressed as 12+log(X/H), where X/H is 
the sum of the ionic abundances. 

a Calculated using the ionic abundances of Ar II and Ar III. 

b Calculatcd using the ionic abundances of Fe II, Fe III, and Fe IV. The ICF applied due 
to the presence of Fe IV is indicated between parentheses. 

c Thc first term gives the measurement error and the error on the ICF if applied. The 
second term gives the error due to the assumed physial conditions f ^4.3^ . 
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Table 8 

Elemental abundances from the optical. 



No/H 


S/H 


Ar/H 


Fe/H 


Reference 


NGC 3603 


8.03 ± 0.11 


7.36 ± 0.08 


6.58 ± 0.17 


6.05 ± 0.10 


(1) 


7.94 ± 0.26 


6.83 ± 0.24 


6.55 ± 0.13 




(2) 


8.08 


7.12 ± 0.09 






(3) 


Average Nc/S 


= 0.82±0.21 








30 Dor 




6.70 ± 0.14 






(3) 


7.83 ± 0.06 


6.99 ± 0.10 


6.26 ± 0.10 


6.39 ± 0.20 


(4) 


7.73 


6.84 


6.14 


6.12 


(5) 


7.88 


6.84 


6.15 ± 0.07 




(6, pos#l) 


7.84 ± 0.08 


6.80 ± 0.09 


6.08 ± 0.05 




(6, pos#2) 


7.84 ± 0.16 


6.82 ± 0.04 


6.21 ± 0.05 




(6, pos#3) 


7.94 ± 0.07 


6.88 ± 0.13 


6.11 ± 0.10 




(6, pos#4) 


7.55 


6.71 


6.03 




(7) 


7.78 ± 0.06 


6.67 ± 0.01 


6.25 




(8) 


Average Nc/S 


= 1.03±0.12 






N66 


7.47 ± 0.01 


6.32 ± 0.02 


5.81 




(6) 


7.24 








(9, N66A) 


7.20 








(9, N66NW) 


7.22 








(10) 


7.22 


6.50 


5.74 




(12, t 2 =0.0013) 


Average Ne/S 


= 1.15±0.02 









Note. — The average Ne/S values were calculated from measurements with 
quoted errors only. 

References. — (1) Garci'a-Rojas ct al. (2006); (2) Tapia ct al. (2001); (3) 
Simpson et al. (1995); (4) Peimbert 2003; (5) Tsamis & Pequignot 2005; (6) 
Vermeij & van dcr Hulst (2002): (7) Rosa & Mathis 1987; (8) Mathis et al. 
1985; (10) Dufour, Shields & Talbot (1982); (11) Dufour & Harlow (1977); 
(12) Peimbert et al. (2000). 



20 



our abundance measurements are less affected by 
electronic temperature determinations, this sug- 
gests that the discrepancy often observed between 
optical recombination lines and CELs is not prin- 
cipally due to temperature fluctuations (see also 
Liu ct al. 2001; Bernard-Salas et al. 2001). 

6.1. MIR abundances 

6.1.1. Neon 

The abundance of neon does not suffer from any 
significant uncertainty due to ionization correc- 
tions f H5.2p . Neon stands out as being the element 
showing the best agreement with the correspond- 
ing optical references ($5). The average abun- 
dance we derive in NGC 3603 (« 7.94±0.11±0.08) 
is relatively close to the value Ne/H=8.08 found by 
Simpson et al. (1995) using the same MIR CELs as 
we observe, but combining different observations 
with possible aperture effects. 

Neon is also the element which shows the small- 
est scatter in its abundance across the various po- 
sitions in each giant H n region. As an illustration, 
the 15 positions in 30Doradus show only 0.11 dex 
dispersion, the 3 positions in NGC 3603 give re- 
markably equal values, while the 11 positions in 
N66 show 0.20 dex dispersion. These findings, to- 
gether with the fact that neon is not expected 
to be depleted on dust grains or incorporated in 
molecules makes it a reliable metallicity tracer, 
contrary to sulfur, argon, and iron (see also £16 . 3[) . 
For this reason, we consider it as a reference for 
our discussion. 

6.1.2. Sulfur 

The average abundances of sulfur we deter- 
mined in the giant H n regions corroborate the 
optical values (Fig. [5} . There might be an indica- 
tion however that our determinations agree best 
with the lowest optically-derived abundances (see 
discussion in §6.2.ip . 

The sulfur abundance shows a larger dispersion 
across each giant H n region than neon. Part of 
the larger dispersion could be attributed a priori 
to ionization corrections due to the presence of S II. 
However, an ICF would change S/H by less than 
0.05 dex (assuming S n/S 111=0.1; §5.1| . which is 
smaller than the total uncertainty on the sulfur 
abundance determination (±0.11 ± 0.08 dex). 



In order to study in more detail the sulfur 
abundance in the H n regions, we compare it to 
the neon abundance (see §6.1.ip . As discussed in 
Thuan et al. (1995), Ne, S, and Ar are products of 
a-processes during both hydrostatic and explosive 
nucleosynthesis in massive stars. These elements 
are thus thought to follow a parallel chemical evo- 
lution in stars except in the case of an extreme 
initial mass function (IMF) because the stars syn- 
thetizing Ne have somewhat larger masses than 
the stars producing S and Ar (Woosley & Weaver 
1995). We observe that Ne/H and S/H trace each 
other remarkably well in the three objects of our 
sample (Fig. [9}. This contrasts with the lack of 
correlation found by Verma et al. (2003) in star- 
burst regions using ISO data, but this is consistent 
with Ne and S production in massive stars. 

The Ne/S ratio appears to be significantly 
larger than the solar value by ~0. 2-0. 3 dex (Fig. [9}. 
Correcting for the presence of S n cannot explain 
this discrepancy f i|5.ip . Note that an extreme IMF 
is also unlikely to explain the sulfur depletion be- 
cause S/H and Ar/H show somewhat different 
behavior with respect to Ne/H although S and Ar 
are produced in similar mass stars f H6.1.3p . 

Interestingly, the Ne/S ratio was found to be 
on average larger in the PNe of the Magellanic 
Clouds than in the Milky Way, which was at- 
tributed to sulfur depletion on molecules and dust 
grains (Marigo et al. 2003; Henry et al. 2004; Pot- 
tasch Sz Bernard-Salas 2006; Bernard-Salas et al. 
2007). This underabundance of sulfur with respect 
to other a-elements has also been observed in H II 
regions in the Milky Way (Simpson et al. 2004) 
and in M83 and M33 (Rubin et al. 2007) with the 
same interpretation. Unlike Si or Fe for instance, 
sulfur is not depleted toward cool diffuse clouds 
(Savage & Sembach 1996). However, in dense re- 
gions and in PDRs, it becomes significantly de- 
pleted (Simpson & Rubin 1990; Verma et al. 2003; 
Pottasch & Bernard-Salas 2006). Goicoechea et al. 
(2006) found that sulfur depletion could be up to 
a factor 4, i.e., 0.6 dex in the Horsehead PDR. For 
even denser regions, it has been established that 
sulfur is about two orders of magnitude more de- 
pleted than C, N, and O in molecular clumps with 
densities n(H)~ 10 3 ~ 5 c m - 3 (Ruffle et al. 1999). 

The position disagreeing the most with the 
Ne/S solar value is 30DOR#8, which shows promi- 
nent silicate absorption, where the ISM is denser 
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Fig. 9.— Top - S/H is plotted against Ne/H. 
'BCDs' stands for blue compact dwarfs, while 
'GHIIR' stands for giant H n regions. Bottom - 
The ratio Ne/S is plotted as a function of Ne/H. 
The gray stripe illustrates the solar proportion be- 
tween Ne and S, based on the range of the latest 
solar abundance determinations. The black circle 
shows the solar abundance ratio. 

and where refractory elements are more likely to 
be depleted onto dust grains. However, the ion- 
ized gas we measure the abundances from is not 
associated with the dense region, as indicated by 
the negligible amount of MIR extinction toward 
30DOR#8 (Lebouteiller et al. in preparation). 
Furthermore, regions showing ratios closer to the 
solar Ne/S (e.g., N66#7, NGC3603#3) show PAH 
emission, which is a sign of relatively dense re- 
gions. Hence, though sulfur depletion onto dust 
grains is likely to occur, there is no trend with 
the spectral charateristics (PAH-dominated, ion- 



ized gas, etc.). This is due to the fact that the 
ionized gas we probe is not entirely associated with 
the PDRs or the embedded regions, because of 
the complex geometry and structure of the lines 
of sight. 

Results in blue compact dwarfs (Izotov & 
Thuan 1999; Wu et al. 2007), and Galactic H n re- 
gions (Simpson et al. 2004) are added for compar- 
ison in Figure [5] BCDs have optical spectra very 
similar to Galactic H II regions, and they sample 
the low-metallicity zone in the plot. The giant H n 
region N 66 in the SMC overlaps the most metal- 
rich BCDs. Note that abundance measurements 
from the optical and MIR studies of BCDs dis- 
agree with each other (crosses and filled circles in 
Figure [9]). However in the case of the BCDs, the 
difference comes from a higher neon abundance 
measured in the MIR, which could be due to the 
fact that the gas probed in the MIR is enriched 
in heavy elements and that sulfur is depleted on 
dust grains. It seems that the trend formed by 
the BCDs and the H n regions of our sample is 
significantly off the solar proportion for metallici- 
ties close to solar, while it barely follows the solar 
proportion at metalllicities lower or equal to that 
of N66. Globally, there might be a hint that the 
Ne/S ratio increases together with the metallic- 
ity, which would be an important proof of sulfur 
depletion on dust grains. 

6.1.3. Argon 

In NGC 3603, Ar/H it is somewhat smaller than 
the optical determination (Fig. [8|). The argon 
abundances in N 66 and 30 Dor are surprisingly 
large even when no ICF accounting for Ar iv is 
applied. The argon abundances determined us- 
ing only the Ar II and Ar in contributions match 
better the solar and BCD values than with ICF 
correction due to the possible presence of Ar IV 
(Fig. [TO]). As an illustration, when using an an ICF 
to account for Ar iv, we find Ar/H in 30 Dor equal 
to the highest solar determinations (while 30 Dor 
has clearly a subsolar metallicity with 0.60 Z Q ; 
^6.31) . These results are consistent with the fact 
that the main ionization stage in the dense gas 
probed in the MIR is Ar III, with a small contri- 
bution due to Ar II ( flOf . 

Argon is produced by oxygen-burning in mas- 
sive stars, and it is expected to follow both sulfur 
and neon nucleosynthesis. Indeed the argon abun- 
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Fig. 10. — Ar/H is plotted against Ne/H. top - 
Ar III ion contribution is solely used, without ion- 
ization corrections. Bottom - An ionization cor- 
rection was applied, due to the presence of Ar iv 
(see q5.3p . See Fig. [5] for a plot description. 

dance Ar/H correlates with Ne/H in the giant H II 
regions (Fig. [TOf , but with a somewhat larger dis- 
persion than the Ne/H vs. S/H correlation. This 
is mostly due to large measurement uncertainties 
on the [Ar n] and [Ar in] line fluxes ( tJ3.5[) . and 
not to the contribution from Ar iv f ^5.3[) . There 
is no position in any of the H n regions showing 
an argon deficiency with respect to neon, even at 
high mctallicitics, which is consistent with argon 
being undepleted on dust grains or in molecules, 
even in the densest regions probed. 

The MIR neon and argon abundance determi- 
nations in NGC 3603 and 30 Dor agree well with 
the solar value (Fig. [P0|l . In contrast, about half 



of the positions in N 66 shows a larger argon abun- 
dance than expected from the solar Ne/Ar. Verma 
et al. (2003) find that argon is overabundant rel- 
ative to neon in a wide variety of star-forming 
galaxies, with metallicities from sub-solar (BCDs) 
to super-solar (WR galaxies). There is no trend 
among the galaxies in their sample. We could see 
the effect of enhanced argon abundance in N66 
but it is more likely due to uncertaities in the 
abundance determination. 

6.1. 4- Iron 

The average Fe/H determinations in NGC 3603 
and 30 Dor agree with each other (Table [7]), 
and barely match Fe/H in the most metal-rich 
BCDs (Fig. [TT]) . The average iron abundance in 
NGC 3603 (6.09 ± 0.37 ± 0.08) compares well with 
the value derived from the optical (6.05 ± 0.10; 
Garcia- Roj as et al. 2006), and implies no differ- 
ential depletion between the gases probed in the 
MIR and in the optical f fl6~2~Tj) . 

The iron abundance in 30 Dor shows a large 
dispersion including the optical value within the 
range. The Fe/H value toward 30DOR#17 is more 
than 3er larger than the average iron abundance 
in 30 Dor. The spectrum of 30DOR#17 shows 
deep silicate absorption probably originating from 
dust associated with an ultracompact H n region. 
In the same spectrum, we also observe the high- 
excitation [O iv] line at 25.89 /im, which is known 
to originate around hot Wolf-Rayet stars (Schacrcr 
& Stasinska 1999) and in shock-heated gas (Lutz 
et al. 1998). Because the region is dominated 
by schocks from the nearby supernova remnant 
(SNR) 30 Dor B (Chu et al. 1992), it is likely that 
the [O iv] line traces gas shocked by the SN. Hence 
the enhanced iron abundance toward 30DOR#17 
is consistent with removal of iron atoms from dust 
grains due to shocks. We plan to investigate in 
more detail the ISM properties toward the SNR 
using the new data provided by an accepted IRS 
Cycle 5 GTO program. 

There are only two Fe/H determinations in N 66 
(positions #1 and #8), both of which are signif- 
icantly smaller than the Fe/H values in the two 
other giant H n regions, and also smaller than 
the BCDs having similar metallicities (as traced 
by Ne/H). This might be due to underestimated 
ionization corrections due to the presence of Fe IV 
in this object ft j5.4[) . However, although N66 is 
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Fig. 11— Fe/H is plotted against Ne/H. See 
Fig. [Ufor a plot description. 

characterized by a globally harder ISRF as com- 
pared to the two other H n regions in our sample, 
the positions N66#l and N66#8 do not show sig- 
nificantly harder ISRFs than all the positions in 
30 Dor. Hence the iron underabundance in N66 
could be genuine. 

The Fe/Ne ratio in the three giant H n re- 
gions is significantly lower than the solar value 
(Fig. [TT|) . This can be explained by iron depletion 
on the surface of very small grains, which are the 
dominant dust component in the ionized gas (see 
e.g., Cesarsky et al. 1996; Verstaete et al. 1996; 
Lebouteiller et al. 2007). Rodriguez (2002) ob- 
served that the [O /Fe] ratio increases with metal- 
licity in BCDs, and explained it by iron depletion 
on dust at high metallicity. We find that the trend 
formed by the BCDs and the giant H n regions 
flattens for large neon abundances, suggesting that 
the iron becomes more and more underabundant 
as the metallicity increases, being consistent with 
iron depletion on dust grains. 

6.2. Dense gas properties 

6.2.1. Depletion onto dust grains 

MIR spectra allow us to probe regions that are 
obscured at optical wavelengths, possibly with dif- 
ferent physical and chemical properties. The av- 
erage sulfur abundances we calculated from IRS 
spectral lines in the giant H n regions corroborate 
the optical determinations (Fig. [S]), even though 



they likely fall towards their lower side. Interest- 
ingly, Simpson et al. (1995) measured the MIR 
[S III] line flux, and used the [S iv] line flux in 
the literature to infer S/H=7.12 in NGC 3603, i.e., 
close to our average determination for this ob- 
ject (6.99 ± 0.11 ± 0.08; Table H). These two 
MIR determinations are significantly lower than 
7.36±0.08, which was obtained from optical lines 
by Garci'a-Rojas et al. (2006). If the relative 
sulfur underabundance in the MIR is real, this 
could hint at a differential sulfur depletion be- 
tween the gases probed in the MIR and in the 
optical. However, it is likely that the underabun- 
dance is driven by uncertainties in the abundance 
determination methods. Ionization correction in 
the MIR could represent 0.04-0.05 dex (gSHJ). On 
the other hand, optical studies do not directly ob- 
serve S iv, which in our case contribute around 
5% to the total sulfur in NGC 3603, and as much 
as 24% in N66 and 17% in 30 Dor. The ICF ap- 
plied in the optical could overcorrect the total sul- 
fur abundance. Morevover, electronic temperature 
variations/uncertainties affect the optical determi- 
nations and to a less extent the MIR determina- 
tions (S3]). 

The study of the refractory element iron sup- 
ports the lack of differential depletion between the 
gases probed in the MIR and in the optical. We 
found that the iron abundance in NGC 3603 and 
30 Dor is similar when derived in the optical and 
in the MIR ( ^6. 1.4ft . This confirms that the pos- 
sible discrepancy seen in S/H between MIR and 
optical observations is not due to additional de- 
pletion in the gas observed in the MIR, but rather 
to uncertainties on the abundance determination. 
This stresses the importance of deriving the abun- 
dance of a refractory element such as iron in dense 
regions. 

6.2.2. Ionization degree 

We found that the chemical composition de- 
rived in the optical and in the MIR is overall 
fairly similar (Fig. [5]). For this reason, it is pos- 
sible to compare the ionic abundances from op- 
tical and MIR wavelengths and interpret possi- 
ble differences in physical conditions. In order to 
trace the gas excitation, we use the ionic abun- 
dance ratio (Ne m/H)/(S in/H). This ratio is 
a good approximation for (Ne m/H)/(Ne n/H) 
in H II regions (Lebouteiller et al. in prepara- 
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tion) and it can be derived in both optical and 
MIR ranges. Based on the values in Table EJ 
the (Ne m/H)/(S m/H) ratio is 3.7 ± 0.7 in 
NGC3603, i.e., about two times smaller than the 
value in the optical (7.7 ± 0.4; Garcia-Rojas et al. 
2006). We find (Ne m/H)/(S m/H) = 7.6 ± 1.5 
in 30 Dor, somewhat smaller than 10.7 ± 0.1 in 
the optical (Pcimbcrt 2003). Finally, we find 
(Ne m/H)/(S m/H) = 9.8 ± 2.0 in N66, as com- 
pared to 10.3 in the optical (Vermeij & van der 
Hulst 2002). Hence the MIR spectra seem to probe 
ionized gas with a degree of ionization equal to or 
lower than that of the gas probed in the optical. 
This could be due to the fact that the SH and LH 
apertures from the IRS do not sample the same 
spatial regions as those observed within narrow 
slits used in the optical. It is also possible that the 
gases probed in the MIR and in the optical have 
different physical conditions (such as density). 

6.3. Metal dispersion and mixing 

6.3.1. Metal enrichment history 

Stellar winds and supernovas explosions from 
the current star-formation episode release newly 
produced elements into the surrounding ISM on a 
short timescale (few 10 6 yr). It is thus natural to 
question the origin of the metals we observe in the 
ionized gas of the giant H n regions. We investi- 
gate in this section the total metal content of each 
giant H n region by comparing their metallicities 
with the Sun and other objects from the same host 
galaxies (H II regions, planetary nebulae). 

First it must be stressed that the solar abun- 
dance determinations have shown strong varia- 
tions over the past few years (see e.g., Pottasch 
& Bernard-Salas 2006). The gray stripes in Fig. [8] 
illustrate the range of solar abundance determi- 
nations for each element over the 1998-2007 pe- 
riod. Neon and argon abundances are particularly 
poorly determined because of the absence of suit- 
able lines in the solar photosphere. Abundance 
determinations of those elements are indirect and 
make use of coronal lines together with a correc- 
tion using other a-elements such as magnesium or 
oxygen. Considering the most extreme determi- 
nations among the recent studies (Asplund et al. 
2006; Feldman & Widing 2003), solar Ne/H de- 
terminations range from 7.84 ±0.06 to 8.08 ±0.06 
while Ar/H ranges from 6.18 ± 0.06 to 6.62 ± 0.06. 



The sulfur solar abundance is better determined, 
with values from 7.14 (Asplund et al. 2006) to 7.33 
(Grevesse et al. 1998). 

We consider neon as being the most reliable 
metallicity tracer available in our study, as com- 
pared to sulfur, argon, and iron f £|6.1.ip . The 
average neon abundance is 7.94(±0.11 ± 0.08) 
in NGC3603, 7.76(±0.11 ± 0.08) in 30 Dor, and 
7.34(±0.11±0.08) inN66 (Table©. Considering 
an average solar value of 7.98, the average neon 
abundance in NGC 3603 implies essentially a solar 
metallicity («0.91 Z ), while 30 Dor is «0.60 Z Q , 
and N 66 is «0.23 Z s . 

The metallicity of NGC 3603 agrees with the 
solar metallicity within uncertainties. This is con- 
sistent with the abundance gradient in the Milky 
Way (see, e.g., Martin-Hernandez et al. 2002) 
given the fact that the galactocentric distance of 
NGC 3603 is almost that of the Sun (8.5 kpc as 
compared to 8 kpc; i j2. 1 .1(1 . It is instructive to 
compare our abundances in 30 Dor and N 66 with 
those in planetary nebulas (PNe) and in other 
H II regions in the LMC and SMC from Bernard- 
Salas et al. (2007). The authors measured abun- 
dances using MIR lines in a similar way to the 
present study so that the comparison does not 
suffer from significant systematic uncertainties. 
Bernard-Salas et al. (2007) found that neon and 
sulfur abundances in PNe agree relatively well 
with those in the H II regions of the LMC and 
SMC, implying that H n regions have not been 
significantly enriched over the past few Gyrs (age 
of the PNe progenitor). The global neon abun- 
dance we inferred from the MIR lines in 30 Dor is 
7.76(±0.11 ± 0.08), remarkably close to the aver- 
age value in the PNe of the LMC <Ne/H>=7.78 
(no quoted errors). The global neon abundance 
we derive in N66 (7.34 ± 0.11 ± 0.08) is consistent 
within error bars with the average <Ne/H> in the 
PNe of the SMC (7.43). In addition, it must be 
noticed that the metallicity in O dwarf stars in 
NGC 346, the stellar cluster associated with N66, 
is 0.2 Z (Bouret et al. 2003), which agrees well 
with our metallicity determination in the ISM of 
the H II region («0.23 Z Q ). 

It seems that the chemical composition of giant 
H II regions compares well with that of PNe and 
young stars, suggesting that the ionized gas of the 
H II regions has not been enriched significantly by 
the current star-formation episode, or by any pre- 
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vious episode more recent than <1 Gyr. If enrich- 
ment has really occured, the metallicity of the H II 
regions has been enhanced by much less than a fac- 
tor 2. Full mixing requires that elements acquire 
the same physical conditions as the ionized gas, 
namely temperature and density, but also viscosity 
and molecular diffusion (Scalo & Elmegreen 2004). 
The effective mixing provided by turbulent diffu- 
sion in H II regions allows metals to mix on spatial 
scales of a few hundreds parsecs over 100 Myr (Roy 
& Kunth 1995; Avillez & Mac Low 2002). Hence 
abundance discontinuities between H II regions in 
a single galaxy should be ubiquitous, especially in 
low-mass galaxies where rotational shear is weak 
(Roy & Kunth 1995). However observations chal- 
lenge this hypothesis and suggest that mixing can 
occur on spatial-scales as large as the galaxy size 
and over time-scales larger than the age of the 
H II regions (see e.g., Kobulnicky 1998; Skillman 
& Kennicut 1993; Noeske et al. 2000; Russel & Do- 
pita 1990). Our results support such spatial- and 
time-scales and question the fate of metals released 
by massive stars in giant H n regions. Following 
the scenario of Tenorio-Tagle (1996), the metals 
could undergo a long cycle in a hot phase before 
mixing with the surrounding ISM. The composite 
IRAC and Chandra image of 30 Dor shows large 
cavities filled with hot gas, surrounded by colder 
and denser ionized gas shells and filaments (Brandl 
et al. in preparation). These are supernova-driven 
outflows generated by the superstellar cluster (see 
also e.g., Redman et al. 2003). The temperature 
of metals in the hot cavities is too high for recom- 
bination to occur, and the density too small for 
collisionally-excited levels to be populated so that 
the abundances we derive from the MIR spectra 
do not account for these newly released elements. 
The superbubbles of hot gas in 30 Dor should re- 
main inside the galaxy for hundreds of Myr accord- 
ing to the scenario of Recchi et al. (2001). It could 
eventually lead to a superwind that will enable the 
escape of metals in the halo (Tenorio-Tagle 1996). 

6.3.2. Small-scale mixing 

Our observations probe the chemical abun- 
dances toward several positions within each giant 
H II region. Neon and sulfur abundances show 
remarkably little dispersion in the three giant H II 
regions (Fig. HJ. Vermeij & van der Hulst (2002) 
analyzed 4 positions in 30 Dor and found neon 



abundances ranging between 7.84 and 7.94. Ana- 
lyzing 15 positions, we find a comparable disper- 
sion, 0.11 dex from the mininum to the maximum 
Ne/H values. In the same object, Vermeij & van 
der Hulst (2002) find sulfur abundances ranging 
from 6.63 to 6.86, and argon abundances rang- 
ing from 6.18 to 6.63, which is a smaller disper- 
sion than we measured (but with fewer positions). 
Our results in NGC 3603 and N 66 also imply that 
the neon abundance and, to a lesser extent, the 
sulfur abundance show little dispersion in these 
regions. Part, if not all, of the dispersion of the 
elemental abundances is due to uncertainties in 
the abundance determinations, such as ionization 
corrections (fjSJ) and the assumed nebula physical 
conditions ( ^4.3|) . In this section, we investigate 
whether there is evidence for small-scale mixing. 

The apparent homogeneity of the abundances 
contrasts with the wide variety of physical re- 
gions each line of sight samples toward each ob- 
ject (PDR, stellar cluster, ionized gas, embedded 
regions). This suggests that the ionized gas from 
which we measure the chemical abundances is not 
necessarily associated with these physical regions. 
In fact, the lack of MIR line extinction toward 
PDRs and toward embedded objects implies that 
the whole nebula is filled with ionized foreground 
gas (Lebouteiller et al. in preparation). Can abun- 
dance inhomogeneities still be inferred? The ma- 
terial ejected from supernovae could reach a hot 
coronal phase before falling back onto the galactic 
disk in the form of molecular droplets (Tenorio- 
Tagle 1996). This material will then be photodis- 
sociated by massive stars and mix with the sur- 
rounding ionized gas. According to this model, 
mixing could occur at very small spatial scales 
(<1 pc). Tsamis & Pequignot (2005) introduced in 
their model of 30 Dor small-scale chemical inhomo- 
geneities at the subparsec-sized scale. The typical 
distance between the IRS observations is — 4pc in 
NGC 3603, ~20pc in 30 Dor, and -15 pc in N66. 
Hence our MIR observations do not have the nec- 
essary spatial resolution to probe such chemical 
inhomogeneities. 

On the other hand, Recchi et al. (2001) pro- 
posed that, mostly because of thermal conduction, 
the SNe type II ejecta start to cool down and mix 
with the cold gas within a few 10 6 yr before the 
formation of hot gas outflows. In this particular 
case, small-scale mixing could be observed locally. 
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In order to assert whether small-scale mixing is 
responsible for the abundance variations we ob- 
serve across each giant H II region, we use Ne/H 
in 30 Dor as a reference. The standard deviation 
of Ne/H is 0.02 dex, implying that the abundance 
of neon could vary as much as ~5% across the 
region. Interestingly, 5% of the neon abundance 
currently present in 30 Dor corresponds to the to- 
tal neon enrichment that the star-forming dwarf 
galaxy IZwl8 experienced so far (Ne/H = 6.40; 
Izotov et al. 1999). Kunth & Sargent (1986) pro- 
posed that the metallicity of IZwl8 represents the 
minimum chemical enrichment of an H II region 
from a starburst episode. It is unlikely that the 
small-scale variations in a single giant H II region 
correspond to the yields from a starburst episode, 
especially considering the fact that the metals in 
the ionized gas are not likely to be cotcmporal with 
the star-formation episode ( §6.3|) . Uncertainties 
on the measurements and on the abundance de- 
termination method (±0.11 ± 0.08 dex for Ne/H, 
i.e., as much as 55%) are too large to constrain 
possible abundance fluctuations of less than 5% 
at the metallicity of 30Doradus. 

7. Conclusions 

We analyzed the chemical abundances in the 
ISM of three giant H II regions, NGC 3603 in the 
Milky Way, 30 Dor in the LMC, and N66 in the 
SMC using the MIR lines observed with the IRS 
onboard Spitzer. 

• Our observations probe the ISM toward vari- 
ous physical regions, such as stellar clusters, 
ionized gas, photodissociation regions, and 
deeply embedded MIR bright sources. The 
spectra show the main ionization stages of 
neon, sulfur, and argon in the ionized gas. 
We also detect [Fe n] and [Fe in] lines. 

• Ionic abundances of Ne II, Ne in, S in, S iv, 
Ar II, Ar in, Fe II, and Fe in were derived. 
The internal variation of electron density 
across a region has no impact on the ionic 
abundance determination. On the other 
hand, we find that electron temperature un- 
certainties and/or intrinsic variations could 
be responsible for an error of 20% at most 
on the abundance determinations. Based on 
the (Ne Iii/H)/(S in/H) ionic abundance ra- 
tio, we find that the optical spectra probe a 



gas with a degree of ionization equal to or 
higher than the gas probed in the MIR. 

• Elemental abundances were determined 
from the ionic abundances. No ionization 
corrections were needed, except for iron. We 
find that neon, sulfur, and argon scale with 
each other, which is expected from stellar 
yields. Abundances do not show any depen- 
dence on the physical region (PDR, stellar 
cluster, embedded region, ...). 

• The Ne/S ratio is larger than the solar value, 
and suggests that sulfur could be depleted 
onto dust grains. The sulfur abundance in 
the MIR agrees best with the lowest optical 
determinations, which is likely due to uncer- 
tainties in the abundance determinations. 

• Iron abundance shows a larger uncertainty 
than Ne/H, S/H, and Ar/H. The compar- 
ison of iron and neon abundances hints at 
significant depletion of iron onto dust grains 
at large metallicities. The agreement with 
the optical determination of Fe/H indicates 
however that there is no differential deple- 
tion on dust grains between the gas probed 
in the MIR and in the optical. 

• Fe/H is found to be spectacularly large in 
one position, corresponding to a supernova 
remnant. This strongly suggest that iron 
atoms have been released from dust grains 
due to schocks from the SN. 

• The metallicity of NGC 3603 agrees with the 
Galactic abundance gradient. The metallic- 
ities of 30 Doradus and N 66 agree well with 
those of the PNe in their respective host 
galaxies. These findings suggest that the gi- 
ant H II regions did not experience a signifi- 
cant metal enrichment for at least 1 Gyr. If 
enrichement occurcd, the metallicity was al- 
tered by less than a factor of two. 

• Neon and sulfur abundances show remark- 
ably little dispersion in the three H n regions 
(e.g., 0.11 dex dispersion in 15 positions in 
30 Dor). Small-scale mixing is apparently ef- 
fective, abundance fluctuations are smaller 
than ^55%. However, internal variations of 
the abundances are likely to be on the or- 
der of <5%, and determining their existence 
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would require a significant improvement of 
the data quality and of the method to be 
evidenced. 
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A. Glossary of abbreviations 

BCD : Blue compact dwarf. 
GEL : Collisionally-excited line. 
GTO : Guaranted time observation. 
ICF : Ionization correction factor. 
IMF : Initial mass function. 
IP : Ionization potential. 
ISM : Interstellar medium. 
ISRF : Interstellar radiation field. 
LMC : Large Magellanic Cloud. 
PAH : Polycyclic aromatic hydrocarbon. 
PDR : Photodissociation region. 
PICS : Photoionization cross-section. 
SMC : Small Magellanic Cloud. 
SNR : Supernova remnant. 



B. Spectra 
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Fig. 12.— Spectra of the MIR sources in NGC3603 (Table [J). The spectral window 20-36 in LH is not 
shown in positions #7, #8, and #9 because of saturated lines and continuum. 
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Fig. 13.— Spectra of the MIR sources in NGC3603 (Table HJ. The spectral window 20-36 in LH is not 
shown in positions #7, #8, and #9 because of saturated lines and continuum. 
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Fig. 14.— Spectra of MIR sources in 30 Dor (Table HJ. 
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Fig. 15.— Spectra of MIR sources in 30 Dor (Table HJ. 
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Fig. 16.— Spectra of MIR sources in 30 Dor (Table HJ. 
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Fig. 17. — Spectra of MIR sources in N66 (TableQ]). Posisions #1, #2, and #3 were observed only with the 
high-resolution modules. 
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Fig. 18. — Spectra of MIR sources in N66 (Table [TJ. Posisions #1, #2, and #3 were observed only with the 
high-resolution modules. 
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